Floquet generation of Majorana end modes and topological invariants 
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We show how Majorana end modes can be generated in a one-dimensional system by varying one 
of the parameters in the Hamiltonian periodically in time. The specific model we consider is a wire 
containing spinless electrons with a p-wave superconducting term and a chemical potential; this 
is equivalent to a spin- 1/2 chain with anisotropic XY couplings between nearest neighbors and a 
magnetic field applied in the z direction. We show that varying the chemical potential (or magnetic 
field) periodically in time can produce Majorana modes at the ends of a long chain. We discuss 
two kinds of periodic driving, periodic (^-function kicks and a simple harmonic variation with time. 
We discuss some distinctive features of the end modes such as the inverse participation ratio of 
their wave functions and their Floquet eigenvalues. For the case of periodic 5-function kicks, we use 
the effective Hamiltonian of a system with periodic boundary conditions to define two topological 
invariants. The first invariant is a well-known winding number while the second invariant has not 
appeared in the literature before. We show that the second invariant correctly predicts the numbers 
of end modes with Floquet eigenvalues equal to +1 and — 1. 

PACS numbers: 03.65.Vf, 75.10.Pq 



I. INTRODUCTION 



Topological phases of quantum matter have been ex- 
tensively studied for several yearsi~— . Typically, these 
are phases which have only gapped states in the bulk 
(which therefore do not contribute at low temperatures 
to properties like transport) but have gapless states at 
the boundaries. (For three-, two- and one-dimensional 
systems, the boundaries are given by surfaces, edges and 
end points respectively). Further, the number of species 
of gapless boundary modes is given by a topological in- 
variant whose nature depends on the spatial dimension- 
ality of the system and the symmetries that it possesses, 
like spin rotation symmetry, particle-hole symmetry and 
time-reversal symmetry. The significance of a topologi- 
cal invariant is that it does not change if the system is 
perturbed (say, by impurities) , as long as the bulk states 
remain gapped and the symmetry of the system is not 
changed by the perturbation. Examples of systems with 
topological phases include two- and three-dimensional 
topological insulators, quantum Hall systems, and wires 
with p-wave superconductivity. 

Recently, there has been considerable interest in sys- 
tems in which the Hamiltonian varies with time in a peri- 
odic way which gives rise to some topological features^—. 
Some of these papers have discussed boundary modes 
and topological invariants^— ' 10 ' 17 ~ 20 ' 22 . Recently a pho- 
tonic topological insulator has been demonstrated exper- 
imentally; a two-dimensional lattice of helical waveguides 
has been shown to exhibit topologically protected edge 
states 2 ^. However, the existence of topological invari- 
ants and the relation between them and the number of 
Majorana modes at the boundary seems to be unclear 
if the driving frequency is small**. Further, the Majo- 
rana boundary modes are of two types (corresponding 
to eigenvalues of the Floquet operator being +1 or —1, 



as discussed below) ; it would be interesting to know how 
the numbers of these two types of modes can be obtained 
from a topological invariant. The effect of time- reversal 
symmetry breaking on the boundary modes have also 
not been studied in detail. In this paper, we address all 
these questions for a one-dimensional model where both 
Majorana end modes and topological invariants can be 
numerically studied without great difficulty. 

The plan of this paper is as follows. In Sec. II we 
present our model. This consists of a lattice model 
of spinless electrons with p-wave superconducting pair- 
ing and a chemical potential. By the Jordan- Wigner 
transformation 2 ^, this is equivalent to a spin-1/2 XY 
chain placed in a magnetic field pointing in the z di- 
rection. Wc discuss the energy spectrum and the three 
phases that this model has when the Hamiltonian is time- 
independent. In Sec. Ill, we review two topological in- 
variants which one-dimensional models with and with- 
out time-reversal symmetry have when periodic bound- 
ary conditions are imposed. In Sec. IV, we discuss our 
numerical method of studying the Floquet evolution and 
the modes which appear at the ends of a system in which 
the Hamiltonian is varying with time in a periodic way. 
In Sec. V, we study what happens when one of the 
terms in the Hamiltonian (the chemical potential in the 
electron language or the magnetic field in the spin lan- 
guage) is given a periodic (5-function kick^ 7 -. We study 
the ranges of parameters in which Majorana end modes 
appear at the ends of an open system and various prop- 
erties of these modes such as their number and Floquet 
eigenvalues. We use the Floquet operator for a system 
with periodic boundary conditions to define two topolog- 
ical invariants. The first invariant is a winding number 
which gives the total number of end modes, while the 
second invariant appears to be new and it correctly pre- 
dicts the numbers of end modes with Floquet eigenvalues 



2 



equal to +1 and — 1. We find that end modes can either 
appear or disappear as the driving frequency is varied, 
and our second topological invariant predicts where this 
occurs. The effect of time- reversal symmetry breaking 
on the end modes is studied; we find that that the end 
modes may survive but they are no longer of the Majo- 
rana type. In Sec. VI, we consider the case in which the 
chemical potential varies in time in a simple harmonic 
way. We summarize our main results and point out some 
directions for future work in Sec. VII. 



II. THE MODEL 

We consider a lattice model of spinless electrons 
which have a uniform nearest-neighbor hopping ampli- 
tude and p-wave superconducting pairing between neigh- 
boring sites, and a site-dependent chemical potential. For 
a finite and open wire with N sites, the Hamiltonian takes 
the form 

JV-l 

H = J2 hUtU+l + fUlU) + A(/n/ n+1 + fl +1 ft)\ 
n=l 
N 

-5>(2/t/„-l), (1) 



where 7, A and /x are all real; we may assume that 7 > 
without loss of generality. In this section we will assume 
that n is time- independent. The operators f n in Eq. ([T]) 
satisfy the usual anticommutation relations {/ m , /«} = 
and {fm, fn} = Smn- (We will set both Planck's constant 
h and the lattice spacing equal to 1 in this paper). We 
introduce the Majorana operators 



a2n-\ = fn + fn and 



i(fn-fl), (2) 



for n = 1, 2, • ■ • , N. We can check that these are Hcrmi- 
tian operators satisfying {a mi a n } = 26 mn . In terms of 
these operators, Eq. (TT|) takes the form 



JV-l 



H — i [ Jx0.2n,a2n+l ~ JyO>2n-l&2n+2] 



N 



^(7 -A) and J y = i( 7 + A). 



(3) 



Note that the Hamiltonian is invariant under the parity 
transformation V corresponding to a reflection of the sys- 
tem about its mid-point, i.e., a 2n — > (— Y) n a-iN+i-2n an( i 

0.271+1 0-2N-2n- 

We can map the above system to a spin- 1/2 XY chain 
placed in a magnetic field pointing in the z direction. We 
define the Jordan- Wigncr transformation from N spin- 



1/2's to 2V Majorana operator 



,26 



<22n-l 



0,2r, 




(4) 



where the c° denote the Pauli matrices at site n, and n = 
1, 2, • • • , N . We then find that Eq. ([3]) can be rewritten 
as 



H 



N-l 

E 

n=l 



r t ~ x ^ 



J n v n v 1 
J y°n°n+l\ 



N 
n=l 



(5) 



The system discussed above is time-reversal symmet- 
ric. The time-reversal transformation involves complex 
conjugating all objects, including i — > —i. With the usual 
convention for the Pauli matrices, Eq. (j4]) implies that 



0-2n 



0-2n 



and a 2n +i -> a 2n +i- 



(G) 



Hence Eq. Q is time-reversal symmetric. 

In order to find the energy spectrum of this system 
in the bulk, we consider a chain with periodic bound- 
ary conditions. We define the Fourier transform fk = 



Sn=i fn£ lkn , where the momentum k goes from 



1 

Vn 

to 7r in steps of 2ir/N. Then Eq. ([1])) can be written in 
momentum space as 



H = E (ft f-k ) 

0<fc<ir 

hk = 2(7 cos k — n) s z + 2Asinfe s v 



fk 

J-k 



(7) 



where the s a are Pauli matrices denoting pseudo-spin. 
(There is a small problem in Eq. (0 for k — and it 
since k = —k in those two cases. We will ignore this 
since these are only two points in the space of all real 
numbers from to tt). The dispersion relation follows 
from Eq. (J7J| and is given b y 28 ' 29 



E k = 



= ^4(7 COS fc — /i) 2 + 



4A 2 sin 2 k. 



(8) 



Depending on the values of 7, A and fj., the system 
has three phases where Ek is non-zero for all values of k. 
The phase diagram is shown in Fig. [1] Phase I lies in the 
region A/7 < and — 1 < ju/7 < 1. In this phase, a long 
and open chain has a zero energy Majorana mode at the 
left (right) end in which a m is non-zero only if m is odd 
(even). This can be seen by considering the extreme case 
J x > and J y = /i = in Eq. ([3]). Then that Hamilto- 
nian is independent of a\ at the left end and a 2 N at the 
right end; hence we have zero energy modes correspond- 
ing to these two operators. In the spin-1/2 language of 
Eq. (0 , phase I corresponds to long-range ferromagnetic 
order of a x . Next, phase II lies in the region A/7 > 
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FIG. 1: Phase diagram of the model in Eq. |T]) as a function 
of fi/'y and A/7. Phases I and II are topological while III is 
non-topological. 



and — 1 < /Lt/7 < 1; here a long and open chain has a zero 
energy Majorana mode at the left (right) end in which 
a m is non-zero only if m is even (odd). In the spin-1/2 
language, this phase corresponds to long-range ferromag- 
netic order of a v . Finally, phase III consists of the two 
regions with fx/'y < — 1 and /u/7 > 1. In this phase, 
there are no zero energy Majorana modes at either end 
of an open chain. In the spin language, this is a paramag- 
netic phase with no long-range order. The three phases 
are separated from each other by quantum critical lines 
where the energy Ek vanishes for some values of k. The 
critical lines are given by /i/7 = ±1 for all values of A, 
and —1 < ///7 < 1 for A = 0. We will see in the next 
section that the three phases can be distinguished from 
each other by a topological invariant which is given by a 
winding number. 



III. TOPOLOGICAL INVARIANTS FOR A 
TIME-INDEPENDENT HAMILTONIAN 

In this section, we review the meaning of a topologi- 
cal phase and the topological invariants which exist for a 
one-dimensional system with a time-independent Hamil- 
tonian which may or may not have time-reversal symme- 
try^. This discussion will be useful for Sec. V where 
we will study if similar topological invariants exist for a 
system in which the Hamiltonian varies periodically with 
time. 

We begin by considering a general Hamiltonian which 
is quadratic in terms of Majorana fermions, 



H = 



2N 

£ 

771,71 — 1 



1 M„ 



(9) 



where M is a real antisymmetric matrix; hence iM is 
Hcrmitian. We can show that the non-zero eigenvalues 



of iM come in pairs ±Xj (where Xj > 0), and the cor- 
responding eigenvectors are complex conjugates of each 
other, Xj and x*. This follows because iMxj = XjXj 
implies iMx* = —XjX*. The zero eigenvalues must be 
even in number and their eigenvectors can be chosen to 
be real. This is because iMxj — implies iMx* = 0, 
and we can then choose the eigenvectors to be the real 
combinations Xj + x* and i(xj — x*). 

Given the time-reversal transformation in Eq. (|6]), we 
see that the Hamiltonian in Eq. © will have time- 
reversal symmetry if the matrix elements M mn are zero 
whenever both m and n are even or both arc odd. Fur- 
ther, let us assume that the system is translation invari- 
ant and has periodic boundary conditions so that M mn 
is only a function of m — n modulo 2N. Defining the 
Dirac fermions f n using Eq. @, we then find that the 
Hamiltonian will have the form given in Eq. (JT]), with^ 



a% k s y + a 3ifc s 



(10) 



where a-i/z.k some real and periodic functions of k. 
The corresponding dispersion is then given by Ek = 

^Ja^ + a 3 fe- Although Eq. ([7]) defines a 2 /z^ omv f° r 
< k < 7r, it is convenient to extend these definitions to 
the entire range — it < k < n. Next we map hk to the 
vector Vk = a2,fcy + &3,fc£ in the z — y plane. Let us define 
the angle <fik = tan -1 (03,^/02^) made by the vector 14 
with respect to the I axis. Following Refs. [l8ll30j . we can 
then define a winding number by following the change in 
4>k as we go around the Brillouin zone, i.e., 



W = 



dk d4>k 
2tt ~~dk 



(11) 



This can take any integer value and is a topological in- 
variant, namely, it docs not change under small changes 
in hk unless hk happens to pass through zero for some 
value of k in which case the winding number becomes 
ill-defined; this can only happen if the energy Ek = 
at some value of k which means that the bulk gap is 
zero. In a gapped phase, therefore, Eq. ([TTj) defines a Z- 
valucd topological invariant. We call a phase topological 
if W 7^ 0; such a phase will have W zero energy Majorana 
modes at each end of long chair>22. If W = 0, the phase 
is non-topological and does not have any Majorana end 
modes. 

We can now look at the three phases discussed after 
Eq. ([5]). We discover, by taking appropriate limits (like 
/i«7, A or /* > 7, A) that the winding number takes 
the values —1, +1 and in phases I, II and III respec- 
tively. 

Next, we note that if time- reversal symmetry breaking 
terms were present in the Hamiltonian in terms pro- 
portional to s x and the identity matrix I will appear in 
hk in addition to terms proportional to s y and s z . Then 
as k goes from — 7r to 71", hk will generate a closed curve 
in three or four dimensions instead of only two dimen- 
sions, and it would not be possible to define a winding 
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number as a topological invariant. However, it turns out 
that one can define a Z 2 -valued topological invariant in 
that cas o 29 i 31 . We find that at k = and it, hk only has 
a component along s z ; this is essentially because k = —k 
in those two cases, hence terms proportional to s x , s y 
and I cannot appear in hk- Let us denote ho = gos z and 
h-n = 9kS z . Assuming that we are in a gapped phase, so 
that hk ^ for all values of k, the ^-valued topological 
invariant is defined as v = sg^g^g^) (here sgn denotes 
the signum function). If v = — 1, the phase is topological 
and has one zero energy Majorana mode at each end of 
a long chain, but if v = 1, the phase is non-topological 
and does not have any Majorana end modes. 

Finally, we can ask what would happen if one con- 
sidered a time-reversal symmetric system which is in a 
topological phase with winding number W (hence hav- 
ing W zero energy Majorana modes at each end of a 
long chain), and introduced a weak time-reversal break- 
ing term in the Hamiltonian. Generally, what happens is 
that pairs of end modes move away from zero energy to 
energies ±E; the number of modes which remain at zero 
energy (and hence are Majorana modes) is 1 if W is odd 
and zero if W is even. Thus, the Z-valucd invariant W 
would reduce to the ^-valued invariant v as v = (— l) 1 * . 

IV. FLOQUET EVOLUTION 

In this section, we study a system in which the Hamil- 
tonian varies periodically in time, namely, the matrix 
M in Eq. @ changes with time as M{t) such that 
M(t + T) = M(t), where T denotes the time period. 

We consider the Heisenbcrg operators a n (t). These 
satisfy the equations 

^ = i[H(t),a n (t)]. (12) 

Given that H(t) = iY. mn a m{t)M mn (t)a n {t), we find 
that 

= 4]T M mn (t)a n {t). (13) 

n— 1 

If a denotes the column vector (aia 2 , • • • , o,2n) T and M 
denotes the matrix M mn , we can write the above equa- 
tion in the form da(t)/dt = AM(t)a{t). The solution of 
this is given by 

a(t) = U(t,0) a(0), 
where U(t 2l h) = Te 4 ^ 2 dtM ^\ (14) 

and 7" denotes the time-ordering symbol. The time 
evolution operator U{ti,t 2 ) is a unitary (in fact, real 
and orthogonal) matrix which can be numerically com- 
puted given the form of M{t). It satisfies the properties 
U(t 2 ,h) = U~ 1 {t 1 ,t 2 ) and U(t 3 ,h) = U(t 3 ,t 2 )U(t 2 ,t 1 ). 

If M(t) varies with a time period T, we will call U(T, 0) 
the Floquct operator. The eigenvalues of U (T, 0) are 



given by phases, e l0j , and they come in complex conju- 
gate pairs if e l8i ^ 1. This is because U(T, 0)ipj = e l9i ipj 
implies that U(T, 0)ip* = e~ lB ^*. For eigenvalues e 1 ^ = 
±1 (these eigenvalues may, in principle, appear with no 
degeneracy), the eigenvectors can be chosen to be real; 
one can show this using an argument similar to the one 
given above for zero eigenvalues of the matrix iM . 

In Sees. V and VI, we will consider two kinds of pe- 
riodic driving of the chemical potential fi(t) with a time 
period T, namely, periodic ^-function kicks2£ and a sim- 
ple harmonic variation with time. In each case, we will 
look for eigenvectors of U (T, 0) which are localized near 
the ends of the chain. Before discussing the specific re- 
sults in next two sections, let us describe our method of 
finding Majorana end modes and some of their general 
properties. 

A convenient numerical method for finding eigenvec- 
tors of U (T, 0) which arc localized at the ends is to 
look at the inverse participation ratio (IPR). We as- 
sume that the eigenvectors, denoted as ipj, are normal- 
ized so that Y^m=i l^ji 171 )] 2 = 1 f° r each value of j; here 
to = 1,2, ••■ , 2N labels the components of the eigen- 
vector. We then define the IPR of an eigenvector as 

Ij = Y^n=i IV- 1 j ' ( m )| 4 - 11 V'j is extended equally over 
all sites so that \t/jj(m)\ 2 = l/(2N) for each m, then 
Ij = l/(2N); this will approach zero as TV — > oo. But 
if ipj is localized over a distance £ (which is of the or- 
der of the decay length of the eigenvector and remains 
constant as N — > oo), then we will have \ipj(m)\ 2 ~ l/£ 
in a region of length £ and ~ elsewhere; then we have 
Ij ~ l/£ which will remain finite as N — > oo. If TV is 
sufficiently large, a plot of Ij versus j will be able to dis- 
tinguish between states which are localized (over a length 
scale <C TV) and states which are extended. Once we find 
a state j for which Ij is significantly larger than 1/(2TV) 
(which is the value of the IPR for a completely extended 
state), we look at a plot of the probabilities \ipj(m)\ 2 ver- 
sus to to see whether it is indeed an end state. Finally, 
we check if the form of \ipj(m)\ 2 and the value of IPR 
remain unchanged if TV is increased. 

In both the periodic driving protocols discussed in 
Sees. V and VI, we find, for certain ranges of the pa- 
rameter values, that U (T, 0) has one or more pairs of 
eigenvectors with substantial values of the IPR. For each 
such pair, we find that the corresponding eigenvalues are 
complex conjugates of each other and they are both close 
to 1 (or —1); the two eigenvalues approach 1 (or —1) as we 
increase the system size TV keeping all the other param- 
eters the same. Let us denote the corresponding eigen- 
vectors by ipi(m) and ip 2 (jn), where m = 1, 2, ■ • • , 2TV. 
In the limit that TV — >• oo and the eigenvalues approach 
1 (or —1), any linear combination of ip\ and ip 2 w ih a l so 
be an eigenvector of U (T, 0) with the same eigenvalue. 
In that limit, suppose that we find that the probabili- 
ties of the two orthogonal linear combinations, given by 
\(ipi(m) ± ^(w))! 2 , are peaked close to to — 1 and 2TV, 
and that they decay as to moves away from 1 or 2TV. 
We can then interpret these linear combinations as edge 
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states produced by the time-dependent chemical poten- 
tial. The deviation of the two eigenvalues from 1 (or —1) 
is a measure of the tunneling between the two edge states. 
The larger the tunneling, the greater is the deviation of 
the eigenvalues from ±1; this, in turn, implies that the 
two edge states decay less rapidly as we go away from 
the ends of the chain since a slower decay increases the 
tunneling between the two states. 

The situation discussed in the previous paragraph 
is similar in some respects to the problem of a time- 
independent double well potential in one-dimension 
which is reflection symmetric about one point, say, x = 0. 
Then the eigenstates of the Hamiltonian are simultane- 
ously eigenstates of the parity operator. The lowest en- 
ergy states in the parity even and parity odd sectors differ 
in energy by an amount which depends on the tunnel- 
ing amplitude between the two wells; the corresponding 
wave functions, denoted by and if)—, are symmetric 
and antisymmetric combinations of wave functions which 
arc localized in the two wells separately. In the limit 
that the tunneling amplitude goes to zero, the two states 
become degenerate in energy; further, the linear combi- 
nations "0+ ± ip— describe states which are localized in 
the two separate wells. In our Floquet problem, the two 
end states with opposite parity have complex conjugate 
eigenvalues of U (T, 0) given by e ±t8 . In the limit that 
TV — > oo and the tunneling between the two states goes 
to zero, the eigenvalues of C/(T, 0) must become degener- 
ate; this can only happen if e ±l6 approach either +1 or 
-1. 

Finally, after finding the end modes, we check if their 
wave functions are real in the limit of large N. We call 
the end modes Majorana if they satisfy three properties: 
their Floquet eigenvalues must be equal to ±1, they must 
be separated by a finite gap from all the other eigenval- 
ues, and their wave functions must be real. 



V. PERIODIC 5-FUNCTION KICKS 

In this section, we consider the case where the chemi- 
cal potential is given (^-function kicks periodically in time. 
One reason for choosing to consider periodic kicks is that 
this is known to produce interesting effects in quantum 
systems such as dynamical localization^!. We will also 
see that this system is considerably easier to study both 
numerically and analytically than the case of a simple 
harmonic timc-depcndcncc which will be discussed in 
Sec. VI. 

We begin by taking the chemical potential in Eq. ^ 
to be of the form 



n(t) 



Co 



ci S(t — nT), 



(15) 



where T — 2tt/lj is the time period and to is the driving 
frequency. Using Eq. ([5]), we note that this system has 
time-reversal symmetry: H*{—t) = H(t) for all values 



of t. (In general, we say that a system has time-reversal 
symmetry if we can find a time to such that H* (to —t) = 
H(t) for all t, and does not have time-reversal symmetry 
if no such to exists). As discussed below, we numerically 
compute the operator U(T, 0) for various values of the 
parameters 7, A, co, ci, u> and the system size TV. Wc 
then find all the eigenvalues and eigenvectors of U (T, 0). 
Since the system is invariant under parity, V ', one can 
choose the eigenvectors of U (T, 0) to also be eigenvectors 
of V. 

The Floquet operator for a periodic <5-function kick can 
be written as a product of two terms: an evolution with a 
constant chemical potential cq for time T followed by an 
evolution with a chemical potential c\8(t — T). Namely, 



U(T,0) 



e 4Mi e 4M T 



(16) 



where M /i are (2iV)-dimcnsional antisymmetric matri- 
ces whose non-zero matrix elements can be found using 
Eqs. © and ®: 



(■Wo)2n+l,2n 
(Mo)2n-l,2n 



— — (Aio)2n,2n+l 

= - (M, 
= - (M, 

= - (M; 



2n+2,2n-l 






CO 


2n,2n-l = 


2 ' 




Cl 


2n,2n-l = 


2 ' 



;(7-A), 
4(7 + A), 



(17) 



for an appropriate range of values of n. However, in order 
to make the time-reversal symmetry more transparent, it 
turns out to be more convenient to use the symmetrized 
expression 

U(T,0) = e 2Ml e 4MoT e 2M \ (18) 

It is easy to show that the Floquet operators in Eqs. (fT!)]) 
and (|18l) have the same eigenvalues, while their eigen- 
vectors are related by a unitary transformation. We will 
see below that the symmetrized form in Eq. (fT5)) leads to 
some simplifications when we derive an effective Hamil- 
tonian and a topological invariant. 

We now consider a 200-site system (hence with a 400- 
dimensional Hamiltonian) with 7 = 1, A = —1, cq = 2.5, 
ci = 0.2 and u = 12. Fig. H] shows the IPRs of the dif- 
ferent eigenvectors. Two of the IPRs clearly stand out 
with a value of 0.142 each. Wc find that they both have 
Floquet eigenvalue e" 10 = —1, and the value of 9 = tt is 
separated by a gap of 0.148 from all the other eigenval- 
ues. The corresponding eigenvectors are localized at the 
two ends of the system and are real; the corresponding 
probabilities are shown in Fig. [31 The state at the left 
end has non-zero a m only if m is even, while the state 
at the right end has non-zero a m only for m odd. It is 
important to note that the periodic driving has produced 
Majorana end modes even though for the parameter val- 
ues given above, fj,(t) > 7 at all values of t according 
to Eq. (fl"5"j) . i.e., even though the corresponding time- 
independent system lies at all times in phase III (dis- 
cussed in the paragraph after Eq. ([5])) which is a non- 
topological phase. 
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FIG. 2: IPRs of different eigenvectors of the Floquet opera- 
tor for a 200-site system with a periodic (5-function kick with 
7 = 1, A = -1, Co = 2.5, ci = 0.2 and lu = 12. The two 
eigenvectors with the largest IPRs both have an IPR equal to 
0.142 and Floquet eigenvalue equal to — 1. 



reason for these features will become clear below. 
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FIG. 4: Plot of the number of end states versus w for a 200- 
site system, with Floquet eigenvalues +1 (n+, blue squares) 
and —1 (n_, red stars), for a periodic S- function kick with 
7 = 1, A = -1, c = 2.5 and ci = 0.2. 
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FIG. 3: Majorana end states for a 200-site system with a pe- 
riodic (5-function kick with 7 = 1, A = — 1, Co = 2.5, Ci = 0.2 
and u — 12. These states correspond to the two eigenvectors 
with the largest IPRs in Fig. H 

We now vary u> to see how many Majorana end modes 
there are at each end of the system, and, more specifi- 
cally, how many of these modes have Floquet eigenvalues 
equal to ±1. We denote the number of eigenvalues lying 
near +1 and —1 by the integers n+ and n_ respectively. 
Fig. [4] shows a plot of n± versus u> for a 200-site sys- 
tem with 7 = 1, A = -1, c = 2.5 and c\ = 0.2. (We 
have checked that these Floquet eigenvalues are sepa- 
rated from all the other eigenvalues by a gap which re- 
mains finite as N becomes large). We see that although 
the number of end modes is not a monotonic function of 
w, the number generally decreases as to increases, even- 
tually becoming zero for a large enough value of lu. The 



Topological Invariants 

We saw above that there are a number of Majorana end 
modes, which can be further separated into n± depending 
on whether the Floquet eigenvalues e 10 lie near +1 or — 1 . 
Further, the eigenvalues 9 = and tt are separated from 
all the other eigenvalues by a gap which remains finite 
as N — > 00. We then expect the integers n + and n_ to 
be topological invariants, i.e., they will not change under 
small changes in the various parameters of the system. 
The only way in which these integers can change is if the 
eigenvalue gap closes and reopens as we vary the system 
parameters. 

We therefore look for a topological invariant for this 
time-dependent problem^ ' 10 ' 18 i 19 . Interestingly, we will 
discover that we can define a topological invariant in two 
different ways: one is a winding number which only gives 
the total number of Majorana modes at each end of a 
chain, and the other also gives the individual values of 
n+ and n~ which are the numbers of end modes with 
Floquet eigenvalues equal to +1 and —1. 

To define the topological invariants, we consider a sys- 
tem with periodic boundary conditions. Then the system 
is translation invariant and the momentum k is a good 
quantum number; the system decomposes into a sum of 
subsystems labeled by different values of k lying in the 
range [0,7r]. For each value of k, we define a Floquet 
operator Uk{T, 0) which is a 2 x 2 unitary matrix. Using 
Eqs. 0, (IIS]) and ([15]). we find that 

Uk{T 0) = e ic i sZ g- i2T [(7cosfe-c )s i +Asinfcs y ] e icis z /^g\ 

where we take k to lie in the full range — 7r < k < tt. 



7 



Let us assume that 2c\/tt is not equal to an integer 
and A/0. We can then prove an interesting fact about 
Uk(T, 0), namely, that it can be equal to ±1 only if fc = 
or 7r and if T is given by a discrete set of values. First, 
given the above conditions on c\ and A, we can show that 
U k (T, 0) ^ ±1 for any value of k ^ or tt. Next, if fc = 
or 7r, we see from Eq. (fT5)) that Uk(T, 0) 7^ ±7 unless 
2T(cq ± 7) + 2ci = ri7r, i.e., unless uj = 2-7r/T satisfies 



47r(c ± 7) 



(20) 



for some integer value of n. The ± sign in Eq. (j20|) corre- 
sponds to k — 7r and respectively Eq. (|2"U|) holds only 
for a discrete set of values ui. For all other values of u, 
therefore, Uk(T, 0) will not be equal to ±7 for any value 
of k. This also means that for all k, the Floquet eigen- 
values (which are given by the eigenvalues of Uk(T,0)) 
will be separated by a gap from ±1. We are now ready 
to define our topological invariants. 

First topological invariant: Given Eq. (|T9j) . let us 
define an effective Hamiltonian /i e //,fc as 



U k (T,0) 



-,—ih. 



(21) 



The structure of Eq. (jTTJ)) is such that /i e //,fc takes the 
form 



a 2 ,k s y + a 3jfc 



(22) 



as in Eq. (|TU|) . (Indeed, this is the reason we choose the 
Floquet operator of the form given in Eq. (jTHJ) rather 
than in Eq. (JTSJ)). Note that Eqs. JTSJ and (gU) do not 
determine h e ff,k uniquely. To define h e ff,k uniquely, we 
impose the condition that the coefficients in Eq. (f2"2"j) sat- 
isfy < \/ ^2 fc + a a i- < n - (It i s possible to impose this 



2,fc "3,*; 

if Z7 fe (T, 0) ^ ±7; this will be true if w does not satisfy 
Eq. (|2"Uj) ). Given the form in Eq. (f2"2"j) . we can then com- 
pute a winding number W as described in Eq. (jlip 

We note in passing that the condition < a\ k 
7r implies that ft e //,fc can be mapped to a point on the 
surface of a sphere whose polar angles (a, /3) are given by 

z 2 fc + a 3 fc ana - ^ = t an ( a 3,fc/ a 2,&)- As fc goes 
from to 2ir, we obtain a closed curve which does not 
pass through the north and south poles. The integer W 
can then be related to the winding number of this curve 
around either the north pole or the south pole. Note 
that the winding numbers around the north and south 
pole are given by the same integer. 

In Fig. we show the closed curves in the (a.2.fc, 03^) 
plane for four values of uj for a 200-site system with 7=1, 
A = — 1, Co = 2.5 and a periodic 5-function kick with 
Ci = 0.2. For uj = 3, 7, 12 and 17, the winding numbers 
around the origin are given by 2, 2, 1 and respectively. 
These agree exactly with the number of Majorana modes 
at each end of an open chain for those values of uj as 
shown in Fig. [51 

In Fig. [51 we compare the number of Majorana modes 
at each end of a chain and the winding number as a 




FIG. 5: Closed curves in the (02, fc, 13, fc) plane for cj = 3, 7, 12 
and 17, for a 200-site system with 7 = 1, A = —1, Co = 2.5 
and a periodic (5-mnction kick with ci = 0.2. The correspond- 
ing winding numbers around the origin (the point (0, 0) shown 
by a red dot) are given by 2, 2, 1 and respectively. 




2 4 6 



FIG. 6: Comparison of the number of Majorana modes at 
each end of a 200-site system (black solid, y-axis on left) and 
the winding number (magenta dashed, ?/-axis on right) as a 
function of ui from 1 to 18, for 7 = 1, A = —1, co = 2.5 
and a periodic 5-function kick with c\ = 0.2. The inset shows 
a range of uj from 0.2 to 1 where there is a large number of 
Majorana modes. 



function of w, for a 200-site system with 7 = 1, A = — 1, 
Co = 2.5 and c\ = 0.2. In preparing that figure, we have 
considered only those values of uj for which Eq. (|20j) is 
not satisfied. We see that the number of end modes and 
the winding number completely agree in the range 0.2 < 
uj < 18. Hence our first topological invariant correctly 
predicts the total number of Majorana modes for a wide 
range of values of uj. 

In Fig. [5J we have not shown the number of Majorana 
end modes for a; < 0.2. For small uj, we see that the 
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number of end modes increases. (We will make this more 
precise below). However, it becomes more and more dif- 
ficult to identify the end modes as lo becomes small; we 
find that there are a large number of what appear to be 
end modes, but many of them have decay lengths which 
are not much smaller than the system sizes that we have 
considered and their Floquet eigenvalues differ slightly 
from ±1. Thus wc have to go to very large system sizes 
to confirm if all of these are really Majorana end modes, 
i.e., if their Floquet eigenvalues approach ±1 and if these 
are separated from all other eigenvalues by a finite gap 
in the limit of infinite system size. 

Second topological invariant: We observe that 
the momenta k = and ir play a special role since 
Uk(T,0) can be equal to ±7 at only those two values. 
Eq. dp) shows that U (T,0) = e"^ ^ and U n (T,0) = 
e Mr& x s ^ wnere we choose &o/tt m the simplest possible 
way, namely, 

4(c -7) 2ci 
oo — + , 

K = + 4 (23) 

LO IT 

where lo = 2-k/T. We now define a finite line segment, 
called L^, which goes from b to b n in one dimension 
which wc will call the z-axis. 

For lo — > oo, the line L u collapses to a single point 
given by z = 2ci/ir. We have assumed earlier that this 
is not an integer. As lo is decreased, L u will move and 
also increase in size. For our system parameters 7 = 1, 
Co = 2.5 and c\ = 0.2, we find that the right end of 
L u , given by b^ in Eq. (f23|) . crosses the point z = n 
with n = 1 at some value of w. At this point, we see 
from Eq. (|T5|) that the Floquet eigenvalue at k = n is 
equal to e" l7r = — 1. Wc therefore expect that when lo 
decreases a little more and L u includes the point z = 1, 
a Majorana mode will appear at each end of an open 
chain with the Floquet eigenvalue equal to —1. For our 
parameters, we therefore predict, by setting b n = 1, that 
the first Majorana end mode will appear at u ~ 16.04. 
This agrees well with Fig.|5]which shows that a Majorana 
end mode first appears in the range 16 < to < 17 and it 
has a Floquet eigenvalue equal to — 1. As a; is decreased 
further, the right end of L u given by b n crosses the point 
z = n with n = 2 at another value of w; Eq. (jTUJ) then 
shows that the Floquet eigenvalue at k = 7r is equal to 
e mn = 1. As ui is decreased a little more, will include 
the point = 2, and we then expect that Majorana end 
modes will appear with the Floquet eigenvalue equal to 
1. For our parameters, b^ — 2 occurs at lo ~ 7.48. This 
also agrees well with Fig. [6] which shows that a Majorana 
end mode appears in the range 7 < to < 8 with a Floquet 
eigenvalue equal to 1. As cj is decreased further, the left 
end of L u , given by bo in Eq. (|2"3")l . crosses the point z = n 
with n = 1 at some value of lo; Eq. (fl9|) then shows that 
the Floquet eigenvalue at k = is equal to e m7r = — 1. 
As lo is decreased a little more, L u no longer includes the 
point z = 1 and we expect that the Majorana end modes 



with Floquet eigenvalue equal to —1 will disappear. For 
our parameters, b n = 1 occurs at lo = 6.88. We see in 
Fig.[S]that a Majorana end mode with Floquet eigenvalue 
equal to —1 disappears in the range 6 < lo < 7. 

The general pattern is now clear. If Co ± 7 are both 
positive, the left and right ends of the line segment L u 
will both move in the +z direction as lo decreases, i.e., as 
T increases. Then a Majorana end mode with Floquet 
eigenvalue (—1)" will appear whenever the right end of 
L u crosses a point z = n, while an end mode with Floquet 
eigenvalue (—1)™ will disappear whenever the left end of 
L u crosses z = n. These will happen, respectively, when 
b n and 60 in Eq. (|23| become equal to an integer n. 

The above arguments can be rephrased as follows. For 
any value of u>, the number of points z = n (where n is 
an integer) which lie inside the line segment L u is equal 
to the number of Majorana modes at each end of a chain. 
Further, the numbers of points with n odd and even will 
give the numbers of end modes with Floquet eigenvalue 
equal to —1 and 1 respectively. We have numerically 
verified these statements for all the values of u shown in 
Fig. SI In Fig. [71 we show bo and b^ (i.e., the left and 
right ends of LJ) as functions of u> for the parameters 
7 = 1, A = —1, Co = 2.5 and c\ = 0.2. The Majorana 
end modes correspond to the integers lying within the 
shaded region. 

It is clear that the numbers of odd and even integers ly- 
ing inside L w are topological invariants since these num- 
bers do not change for small changes of the system pa- 
rameters. These numbers can change only at values of to 
where either bo or b^ in Eq. (|23p becomes equal to an in- 
teger. When that happens, Eq. (|T9|) becomes equal to ±7 
at either k = or n, and there is no gap to the Floquet 
eigenvalues at neighboring values of k. 

Finally, we note that as u> —> 0, the number of integers 
(and hence the number of Majorana modes at each end of 
a chain) lying between bo and b^ will approach 87/w. Wc 
see in Fig. [51 particularly in the inset, that the number 
of end modes does diverge as 8/lo (recall that we have 
set 7=1). 



Effect of Time-reversal Symmetry Breaking 

Given a periodically driven time-reversal symmetric 
system which has Majorana end modes (namely, modes 
with real eigenvectors and Floquet eigenvalues equal to 
±1 which are separated from all other eigenvalues by a 
gap), we may ask what would happen if we add small 
terms which break time-reversal symmetry. We discover 
that the end modes persist with their Floquet eigenvalues 
continuing to be separated from all other eigenvalues by 
a gap. However, the Floquet eigenvalues move slightly 
away from ±1 in complex conjugate pairs, and the eigen- 
vectors become complex; hence they can no longer be 
called Majorana modes. This is illustrated in Fig.[5]which 
shows the Floquet eigenvalues for the modes at each end 
of the chain for the time-reversal symmetric case given 
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to 

FIG. 7: Plot of bo and as a function of u) for a system with 
7 = 1, A = —1, co = 2.5 and a periodic 5-function kick with 
ci = 0.2. For each value of lo, the number of even and odd 
integers lying in the shaded region between bo and b n gives 
the number of Majorana modes at each end of a chain with 
Floquet eigenvalues equal to +1 and —1 respectively. 

in Eq. (fPo)) . while Fig. [9] shows the Floquet eigenvalues 
for a case with 

oo 

Mn(*) = c + ci ^2 S(t-nT) 

71 — — OO 

+ c 2 J2 S(t---nT), (24) 

n=— oo 

which breaks time-reversal symmetry. Although we can- 
not clearly see from Fig. [5] that the Floquet eigenvalues 
of the end modes have moved away from ±1, we have 
checked numerically that this is so. For a 2800-site sys- 
tem (this is a large enough system size that there is no 
mixing between the two ends), we find that at each end, 
the Floquet eigenvalues near —1 are given by — l±0.0015i 
and — 1 ± 0.0037i, and the eigenvalues near +1 are given 
by 1 ± 0.0007i and 1 ± 0.004CM. 

VI. SIMPLE HARMONIC VARIATION WITH 
TIME 

In this section, we discuss the case where the chemi- 
cal potential is independent of n but varies harmonically 
with t. Namely, the chemical potential in Eq. ((3]) is of 
the form 

li(t) = c + ci cos(wt + </)). (25) 

The Floquet operator can be written as the time-ordered 
product 

U{T,0) =Te 4 fo dtM W ) ( 2g ) 




real part of eigenvalues 



FIG. 8: Floquet eigenvalues close to ±1 for a 1000-site system 
with 7 = 1, A = —1, co = 2.5, and a (5-function kick with ci = 
0.2 at t = which is repeated with a time period T = 2n. For 
this time-reversal symmetric case, there are four eigenvalues 
at exactly +1 and —1 each, separated by a gap from all other 
eigenvalues; these are shown more clearly in the inset. 




real part of eigenvalues 



FIG. 9: Floquet eigenvalues close to ±1 for a 1000-site system 
with 7 = 1, A = —1, co = 2.5, and two (5-function kicks with 
ci = 0.2 and C2 = 0.1 at t = and T/4 which are repeated 
with a time period T = 2n. For this case with no time-reversal 
symmetry, there are four eigenvalues close to but not exactly 
at +1 and —1, separated by a gap from all other eigenvalues; 
these are shown more clearly in the inset. 
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where M(t) is an antisymmetric matrix with the non-zero 
elements 

(M) 2n+h 2n = - (M) 2 „,2n+i = -4(7 -A), 

(A0an-l,2n+a = - (M) 2 „+2,2n-l = - 7(7 + A), 
(-Mj2n-l,2n = _ {M)2n,2n-l 

= -[c + cicos(ut + 4>)). (27) 

Unlike the case of the periodic <5-function kick, the Flo- 
quet operator is no longer a product of only two opera- 
tors; it has to be computed by dividing the time-period 
T into a large number of time steps of size At each, and 
then multiplying T/ At operators together. Finally, we 
have to check that the results do not change significantly 
once At has been made sufficiently small. Hence this 
problem takes much more computational time. For the 
same reason, a numerical calculation of the Floquet oper- 
ator Uk (T, 0) takes more time here than the correspond- 
ing expression given in Eq. (|19[) for a periodic (5-function 
kick; hence the topological invariant is also harder to cal- 
culate. 

Having computed the operator U(T,0), where T = 
27r /u, we again find all the eigenvalues and eigenvectors 
of U (T, 0) which are also eigenvectors of the parity op- 
erator V . As functions of the parameters 7, A, cq, c±, 
<j>, lu and N, we find that the qualitative features of the 
Majorana end modes that we find are similar to the case 
of the periodic (5-function kicks. As before, we find that 
end modes can appear even when the chemical potential 
places the corresponding time-independent system in a 
non-topological phase at all times t. 

Finally, we discuss the effect of the phase (j> in Eqs. (|25|) 
and (|27|) . The Floquet operator U(T, 0), now denoted by 
U^iTyO), clearly depends on <f>. However, we can show 
that the eigenvalues of U$ (T, 0) are independent of (jjM. 
To see this, note that a shift in the phase <p by an amount 
5 is equivalent to a shift in time by the amount 8/ui. 
Hence 

Efy(T,0) = U {T + ^/lo^/lo) 

= U (T + J>/u;,T)Uo(T,ct>/uj) 

= [7 o (T + 0/ W ,T)C/ o (T,O)C/ o - 1 ( ( /)/ W ,O) 

= U {4>/uj,Q)U {T,Q)U^ x {cj>/u),Q), (28) 

where we have used the fact that Uq(T + <fi/uj,T) = 
U (<f>/uJ,0). Eq. (EHJ) shows that [/^(T,0) is re- 
lated to Uq(T, 0) by a unitary transformation involv- 
ing {7o(</>/w, 0); hence they have the same eigenvalues 
while their eigenvectors are related by the same unitary 
transformation. This also implies that studying how an 
eigenvector corresponding to a particular eigenvalue of 
U^,(T, 0) changes with (j> is equivalent to studying how 
that eigenvector changes with time under evolution with 

The effect of a phase change on the wave functions 
of the Majorana end modes can sometimes be quite 



dramatic. We consider a 100-site system with 7=1, 
A = -1, c = 2.5, ci = 0.3 and u = 14. Figs. [TO] and 
ITTI show the probabilities of the two end modes for <f> = 
and 7r/2 respectively; equivalently, we can think of these 
figures as showing the effect of evolving the end modes 
by a time of T/4. We see that the detailed form of the 
end mode wave functions are quite different in the two 
cases. For = 0, the wave function a m of the mode at 
the left (right) end is non-zero only if m is even (odd). 
For <j> = t/2, both end modes have wave functions in 
which a m is non-zero for both even and odd values of m. 
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FIG. 10: Plot of the probabilities of the two end states for a 
100-site system with 7 = 1, A = —1, co = 2.5, and a simple 
harmonic driving with ci = 0.3, ui = 14 and <f> = 0. 
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FIG. 11: Plot of the probabilities of the two end states for a 
100-site system with 7 = 1, A = —1, Co = 2.5, and a simple 
harmonic driving with ci = 0.3, ui = 14 and <j> = tt/2. 
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VII. CONCLUSIONS 

In this work we have shown that periodic driving of 
a one-dimensional model of electrons with p-wave super- 
conductivity or a spin-1/2 XY chain can generate Majo- 
rana modes at the ends if we have a large and open sys- 
tem which is time-reversal symmetric. To simplify the 
calculations, we have mainly studied the case in which 
the chemical potential of the electrons (or the magnetic 
field in the z direction) is given a periodic 5-function kick. 
However, similar results are found when the chemical po- 
tential (or magnetic field) is driven in a simple harmonic 
way. 

The Majorana end modes exist only for very large sys- 
tem sizes and have three characteristic features: the Flo- 
quet eigenvalues are exactly equal to ±1, they are sepa- 
rated from all the other eigenvalues by a finite gap, and 
the wave functions are real. If the system is not time- 
reversal symmetric, we find that there may still be end 
modes whose eigenvalues are separated from all the oth- 
ers by a finite gap may; however, the eigenvalues are no 
longer exactly at ±1, and the wave functions are not real. 
Hence these cannot be called Majorana modes. 

In analogy with the known topological invariants which 
predict the number of zero energy Majorana end modes 
for a system with a time-independent Hamiltonian, we 
have studied if the driven system has topological in- 
variants which can correctly predict the number of end 
modes. We have shown that there are two topological 
invariants which work for a wide range of the driving 
frequency uj for the case of the periodic (5-function kick. 
The first invariant is a winding number which is similar in 
form to the topological invariant for a time-independent 
Hamiltonian with time-reversal symmetry; this invariant 



gives us the total number of end modes. The second in- 
variant is more sensitive in that it can separately give 
us the numbers of end modes with Floquet eigenvalues 
equal to +1 and —1. The second invariant also gives us 
a simple condition which can predict the values of to at 
which end modes appear or disappear. 

Recently there has been considerable excitement over 
claims of the detection of Majorana modes in semicon- 
ducting/superconducting nanowires^— following some 
theoretical proposals^—. A zero bias peak has been 
observed in the tunneling conductance into one end of 
the nanowire, and it has been suggested that this is the 
signature of a Majorana end mode. Our results can 
be tested in similar systems by applying a gate volt- 
age to the nanowire which varies periodically in time 
in some way. One would like to see if such a time- 
dependent gate voltage can give rise to a zero bias peak 
in a system which would be in a non-topological phase 
in the absence of such a voltage. An important question 
which needs to be investigated in this context is how 
the Majorana end modes appear in the steady state af- 
ter the oscillatory part of the gate voltage is switched 
on. This would require a treatment of relaxation mecha- 
nisms such as electron-phonon interactions^. Finally, the 
effects that disorder in the chemical potential^— and 
electron-electron interactions^— have on the Majorana 
end modes also need to be examined. 
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